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Abstract: We study transport phenomena in p-wave superfluids in the context of 
gauge/gravity duahty. Due to the spacetime anisotropy of this system, the tensorial struc- 
ture of the transport coefficients is non-trivial in contrast to the isotropic case. In partic- 
ular, there is an additional shear mode which leads to a non-universal value of the shear 
viscosity even in an Einstein gravity setup. In this paper, we present a complete study of 
the helicity two and helicity one fluctuation modes. In addition to the non-universal shear 
viscosity, we also investigate the thermoelectric effect, i.e. the mixing of electric and heat 
current. Moreover, we also find an additional effect due to the anisotropy, the so-called 
fiexoelectric effect. 
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1. Introduction 



Hydrodynamics is a very powerful description of systems close to equilibrium. Its focus is 
on slowly varying fluctuations with frequency oj and momentum k smaller than the typical 
length scale, the mean free path. Hydrodynamics may be seen as the low-energy effective 
description of interacting systems. Gauge/gravity duality is a very useful tool to further 
develop the hydrodynamic description for various systems. New transport phenomena 
have been uncovered by studying systems which violate parity by an anomaly [1-3]. The 
transport in a system which shows the chiral magnetic effect induced by an axial anomaly 
has been studied in [4-7]. Effects of anisotropy in strongly coupled systems have been 
discussed in [8,9]. Recently the hydrodynamic description for s-wave superfluids which 
may violate parity has been investigated in [10-12]. 

The hydrodynamical description of superfluids is interesting since an Abelian symme- 
try is spontaneously broken. Due to the spontaneous breaking of a continuous symmetry, a 
Nambu-Goldstone boson appears in the spectrum. Since it is massless, it behaves as hydro- 
dynamic mode and has to be included into the hydrodynamical description. In this paper 
we study p-wave superfluids where in addition to the Abelian symmetry, the rotational 
symmetry is spontaneously broken and thus more Nambu-Goldstone bosons appear in the 
spectrum. This leads to an anisotropic fluid in which the transport coefficients depend on 
the direction, i. e. they are tensors. In the case we study here the fluid is transversely sym- 
metric, i. e. the system has an SO{2) symmetry and wc can use this symmetry to reduce 
the tensors to the minimal amount of independent quantities. For instance, the viscosity 
which relates the stress T^*^ in a fluid with the strain V^Up + ^ pU\ given in terms of the 
four velocity of the fluid tt^ is parametrized by a rank four tensor ^q^^^■^P (see e. g. [13, 14]). 
Using the symmetry we find two independent shear viscosities, in contrast to only one in 
the isotropic case, i. e. 50(3) symmetry. 

The shear mode is the transversely polarized fluctuation given for instance by VyUz + 
V zUy for a momentum in x direction. In the isotropic case this is the unique shear mode 
since any momentum can be rotated into the x direction by the 5*0(3) rotational symmetry. 
In the transversely symmetric case, two momenta, one along and one perpendicular to the 
favored direction, e. g. the x direction, must be considered. Thus there are two shear modes. 
If the momentum is along the favored direction, the situation is similar to the isotropic 
case and the strain is again VyUz + V zUy. However if the momentum is perpendicular to 
the favored direction say in y direction, the situation changes dramatically. Now the little 
group is given by the discrete group Z2 and the strain is given by VxUz + ^ zUx- Since the 
shear viscosity can be evaluated at zero momentum, we can characterize the fluctuations 
with respect to the full symmetry group which is in the transversely symmetric case SO{2). 
The first fluctuation VyUz + VzUy is a helicity two state as the shear mode in the isotropic 
case is. The second fluctuation V^Uz + '^zUx however transform as helicity one state under 
the rotational symmetry. This transformation property is due to the rotational symmetry 
breaking and will be very important in this paper. 

In the context of gauge/gravity duality, the spontaneous breaking of continuous sym- 
metries by black holes developing hair was first achieved in [15] and later used to construct 
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holographic superconductors/superfluids by breaking an Abehan symmetry [16,17]. Along 
this line also p-wave superconductors/superfluids have been constructed [18] and gave rise 
to the first string theory embeddings of holographic superconductors/superfluids [19-21]. 
In order to obtain the effects of spontaneous rotational symmetry breaking in the hydro- 
dynamics of p-wave superfluids, we have to take the back-reaction of the superfluid density 
into account, i. e. we consider the effect of the superfluid density on the energy- momentum 
tensor. This was obtained e.g. in [22]. On the gravity side, the p-wave superfluid state 
corresponds to an asymptotically AdS black hole which carries vector hair. 

A very famous result in the context of gauge/gravity duality is that the ratio between 
shear viscosity and the entropy density is universal [23-25]. The ratio is the same for all 
field theories which have a Einstein gravity dual, i. e. the field theory is a large gauge 
theory at infinite 't Hooft coupling A. This universality can be proven as follows: The 
shear mode is the only mode which transforms as a helicity two mode under the SO{2) 
little group and thus decouple from all the other modes. In addition it can be shown that 
the low energy dynamics of this mode is trivial such that the ratio is completely determined 
by gravitational coupling constant which is universal. The universality is lost if finite A^ 
and/or coupling is considered for instance by adding a Gauss-Bonnet term to the gravity 
action (see e.g. [26,27]). 

In the letter [28] we have shown that universality is also absent even at leading order 
in N and A if the fluid is anisotropic. In this case, the universality is lost since one of the 
different shear modes transforms as a helicity one mode under the rotational symmetry 
and can therefore couple to other helicity one modes present in the system. The coupling 
generates non-trivial dynamics which lead to a non-universal behavior of the shear vis- 
cosity. This result is valid for a field theory dual to Einstein gravity without additional 
contributions to the gravity action. In this paper we present the detailed calculations for 
this result. This calculation was suggested already in [29]. We study the complete set 
of the helicity two and one modes in the SU{2) Einstein- Yang-Mills theory in the broken 
phase at zero momentum. 

Along this calculation we find some additional transport phenomena: the thermoelec- 
tric effect in the transversal directions and the fiexoelectric effect. The thermoelectric effect 
is the phenomenon that the electric and heat current mix since charged object transport 
charge as well as energy. This effect has already been studied for holographic s-wave su- 
perfluids [17,30]. We find that the thermoelectric effect in the transversal directions agrees 
with the result found for s-wave superfluids. The flexoelectric effect is known from nematic 
liquids which consists of molecules with non-zero dipole moment (see e. g. [14]). A direction 
can be preferred by the dipoles. In this anisotropic phase, a strain can lead to effective 
polarization of the liquid and an electric field applied to the liquid can lead to a stress. 
This is the first appearance of this effect in the context of gauge/gravity duality. 

The paper is organized as follows: In section 2 we review the holographic setup in 
which p-wave superfluids are constructed and describe their behavior in equilibrium. In 
section 3 we study perturbations about equilibrium. We characterize the fluctuations in 
terms of their transformation under the symmetry groups and determine their equations of 
motion. In addition we calculate the on-shell action and read off the correlation functions. 



-3- 



In section 4 we extract the transport properties out of the correlation functions and find 
the non- universal shear viscosity, the thermoelectric effect and the flexoelectric effect. We 
conclude in section 5. In the appendix A we discuss holographic renormalization. The 
gauge covariant fields are constructed in appendix B. In appendix C we review the numer- 
ical evaluation of correlator when operator mixing is present. Some general remarks on 
anisotropic fluids are given in appendix D. 



2. Holographic Setup and Equilibrium 

We consider SU{2) Einstein- Yang-Mills theory in (4 -|- l)-dimensional asymptotically AdS 
space. The action is 



S= Id^xy/^ 



2Ki ^ ' 4^^ 



+ 'S'bdy , (2.1) 



where is the five-dimensional gravitational constant, A = — -p- is the cosmological con- 
stant, with L being the AdS radius, and g is the Yang-Mills coupling constant. The SU{2) 
field strength F'^,^^ is 

Ft,M = dMA% - d^Alj + e-'-A\,A% , (2.2) 

where capital Latin letter as indices run over {t,x,y, z,r}, with r being the AdS radial 
coordinate, and e"'"^ is the totally antisymmetric tensor with e^^^ = -|-1. The A'^^ are the 
components of the matrix-valued gauge field, A = A'^r'^dx'^ , where the t"" are the SU{2) 
generators, which are related to the Pauli matrices by = cr"/2i. 5bdy includes boundary 
terms that do not affect the equations of motion, namely the Gibbons-Hawking boundary 
term as well as counterterms required for the on-shell action to be finite. Finally it is 
convenient to define 

a = ^ , (2.3) 
9 

which measures the strength of the backreaction. 

The Einstein and Yang- Mills equations derived from the above action are 

4 2 (rr. 1^ P, 



Rmn + J29MN = /^s [Tmn — ^Tp qmn j , (2.4) 

y^jpaMN ^ -e^^cj^^^pcMN ^ (2.5) 

where the Yang- Mills stress-energy tensor Tmn is 

Tmn = ^ (^F^mF^'^n - JsAfivi^^QF'^^^) • (2.6) 
2.1 Hairy Black Hole Solution 

Following ref. [18,22], to construct charged black hole solutions with vector hair we choose 
the gauge field ansatz 

^ = (^(r)r^dt + «;(r)rMx . (2.7) 
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The motivation for this ansatz is as foUows. In the field theory we will introduce a chemical 
potential for the U{1) symmetry generated by r^. We will denote this U{1) as L'^(l)3. The 
bulk operator dual to the U{l)s density is A^, hence we include A^(r) = (f){r) in our 
ansatz. We want to allow for states with a nonzero {Jf), so in addition we introduce 
A]^{r) = w{r). With this ansatz for the gauge field, the Yang-Mills stress-energy tensor in 
eq. (2.6) is diagonal. Solutions with nonzero w{r) will preserve only an SO{2) subgroup of 
the 50(3) rotational symmetry, so our metric ansatz will respect only 50(2). In addition 
the system is invariant under the Z2 parity transformation Py: 2; — )• —x and w — )■ —w. 
Furthermore, given that the Yang-Mills stress-energy tensor is diagonal, a diagonal metric 
is consistent. Our metric ansatz is [22] 

ds^ = -N{r)airfdf + ^T^dr^ + r'^f(r)-^dx^ + r'^firf (d/ + dz^) , (2.8) 

N(r) ' 

withiV(r) = -?^ + ^. For our black hole solutions we will denote the position of the 
horizon as r/j. The AdS boundary will be at r —t- 00. 

Inserting our ansatz into the Einstein and Yang-Mills equations yields five equations of 

motion for m(r), cr(r), /(r), (/>(r), w{r) and one constraint equation from the rr component 

^^ 

dr) 



of the Einstein equations. The dynamical equations can be recast as (prime denotes ^ 



O" = „ ..o h cr — h 



3riW 1/2 3r 

^ _(lfw^ a^fw'^ _.(^_f.N^.^\ (2.9) 
J 3r2Ar2^2 + 3^2 J f ^ N a 



The equations of motion are invariant under four scaling transformations (invariant 
quantities are not shown), 

(/) 0" — )• Act , (/)—)• X(p, 

{II) f^Xf, w^X-^w, 

(HI) r ^ Xr , m — )• X^m , w — Xw , (f> Xcp , 

{IV) r ^ Xr , m — )• X^m , L — XL , — ^ , a — )■ Aa , 

A 

where in each case A is some real positive number. Using (I) and (II) we can set the 
boundary values of both a{r) and /(r) to one, so that the metric will be asymptotically 
AdS. We are free to use (III) to set to be one, but we will retain as a bookkeeping 
device. We will use (IV) to set the AdS radius L to one. 



r^N \r a 

wcjP' , (I 4/' N' 
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A known analytic solution of the equations of motion is an asymptotically AdS Reissner- 
Nordstrom black hole, which has (^(r) = fj, — q/r'^, w{r) = 0, a{r) = f{r) = 1, and 
N{r) = (r^ — + j , where uiq = + and q = jj,rf^. Here n is the value of 



0(r) at the boundary, which is the U{l)s chemical potential in dual field theory. 

To find solutions with nonzero w{r) we resort to numerics. We will solve the equations 
of motion using a shooting method. We will vary the values of functions at the horizon until 
we find solutions with suitable values at the AdS boundary. We thus need the asymptotic 
form of solutions both near the horizon r = Vh and near the boundary r = oo. 

Near the horizon, we define eh = — — 1^1 and then expand every function in powers 
of €h with some constant coefficients. Two of these we can fix as follows. We determine Vh 
by the condition N{rh) = 0, which gives that m{rh) = r\l2. Additionally, we must impose 
■^{^h) = (f'ifh) = for ^ to be well-defined as a one- form (see for example ref. [31]). The 
equations of motion then impose relations among all the coefficients. A straightforward 
exercise shows that only four coefficients are independent, 



where the subscript denotes the order of (so cjg is the value of a{r) at the horizon, etc.). 
All other near-horizon coefficients are determined in terms of these four. 



in powers of with some constant coefficients. The equations of motion again impose 
relations among the coefficients. The independent coefficients are 



where here the subscript denotes the power of e^,. All other near-boundary coefficients are 
determined in terms of these. 

We used scaling symmetries to set cTq = /q = 1. Our solutions will also have = 
since we do not want to source the operator J'^ in the dual field theory (J7(l)3 will be 
spontaneously broken). In our shooting method we choose a value of /x and then vary 
the four independent near-horizon coefficients until we find a solution which produces the 
desired value of /x and has (Jq = /q = 1 and Wq = 0. 

In what follows we will often work with dimensionless coefficients by scaling out factors 
of rh- We thus define the dimensionless functions rh{r) = m{r)/r^, (l){r) = (p{r)/rh and 
w{r) = w{r)/rh, while /(r) and a{r) are already dimensionless. 

2.2 Thermodynamics 

Next we will describe how to extract thermodynamic information from our solutions [22]. 
Our solutions describe thermal equilibrium states in the dual field theory. We will work in 
the grand canonical ensemble, with fixed chemical potential fi. 

We can obtain the temperature and entropy from horizon data. The temperature T is 
given by the Hawking temperature of the black hole. 




(2.10) 



Near the boundary r = oo we define e;, = (^) <C 1 and then expand every function 




(2.11) 




(2.12) 
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is the surface gravity of the black hole, with ^ being the norm of 



rh 



Here k = y/duid^i 

the timelike Killing vector, and in the second equality we write T in terms of near-horizon 
coefficients. In what follows we will often convert from to T simply by inverting the 
above equation. The entropy S is given by the Bekenstein-Hawking entropy of the black 
hole, 



S 



2TrV 



^3 



27r4 



(2.13) 



where Af^ denotes the area of the horizon and V = J d^x. 

The central quantity in the grand canonical ensemble is the grand potential Q. In 
AdS/CFT we identify Q with T times the on-shell bulk action in Euclidean signature. We 
thus analytically continue to Euclidean signature and compactify the time direction with 
period 1/T. We denote the Euclidean bulk action as / and /on-sheii as its on-shell value (and 
similarly for other on-shell quantities). Our solutions will always be static, hence /on-sheii 
will always include an integration over the time direction, producing a factor of 1/T. To 
simplify expressions, we will define / = I/T. Starting now, we will refer to / as the action. 
/ includes a bulk term, a Gibbons-Hawking boundary term, and counterterms. 



I = hulk + IgH + -fcT ■ 



(2.14) 



^bulk''*^^^ and IqiI^'^^^ exhibit divergences, which are canceled by the counterterms in Iqt- 
To regulate these divergencies we introduce a hypersurface r = rbdy with some large but 
finite rbdy We will ultimately remove the regulator by taking rtdy oo. For our ansatz, 
the explicit form of the three terms may be found in [22]. Finally, CI is related to the 
on-shell action, /on-sheih a,s 

^ ~ -^on-shell- 



The chemical potential fj, is simply the boundary value of A^^r) = 
density (J^) of the dual field theory can be extracted from /on-sheii by 



(2.15) 
The charge 



iJi) 



lim 



61 



on-shcU 



3^,2 



27r'^a 



3„h^ 



(2.16) 



Similarly, the current density (i7f ) is 

'^-^on-shcU 



lim 



rbdy-*-oo Ml (rbdy) 



27r3a2_ 



12cj^ 



2 n3 



(2.17) 



The expectation value of the stress-energy tensor of the GET is [32, 33] 



/T- \ 1. 2 (^ion-shell i. 

{%v) = hm — ^ = lim 



=''bdy 



(2.18) 
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where small Greek letter as indices run over the dual field theory directions {t, x, y, z} and 
Kfj_i, = N{r) drjfxu is the extrinsic curvature. We find 

irtt) = 3-2 r^- — rhl , 



CT"-) = -2r'7 ^ K - 8/2^ , (2.19) 



Ua^^ - (0t) a 

'5 (l2a^^ - (0^)^a2 

Notice that (Tt^) = {Tty) = {Ttz) = 0. Even in phases where the current {Jf) is nonzero, 
the fluid will have zero net momentum. Indeed, this result is guaranteed by our ansatz 
for the gauge field which implies a diagonal Yang-Mills stress-energy tensor and a diagonal 
metric. The spacetime is static. 

Tracelessness of the stress-energy tensor (in Lorentzian signature) implies {Ttt) = 
{Txx) + {Tyy) + {Tzz), whlch Is iudccd true for eq. (2.19), so in the dual field theory 
we always have a conformal fluid. The only physical parameter in the dual field theory is 
thus the ratio fi/T. 

For rfiQ = 2 + ^M^' ^0 = 1; '^i = /2 = 0, and (f>\ = —jl we recover the correct 
thermodynamic properties of the Reissner-Nordstrom black hole, which preserves the 50(3) 
rotational symmetry. For example, we find that (Txx) = {Tyy) = {Tzz) and Q = —V{Tyy), 
i.e. {Tyy) is the pressure P. For solutions with nonzero (iTf), the S0{3) is broken to 
50(2). In these cases, we find that (Txx) 7^ {Tyy) = (Tzz)- In the superfluid phase, both 
the nonzero (Jf) and the stress-energy tensor indicate breaking of 50(3). Just using the 
equations above, we also find 

^ = -y>Jyy) ■ (2-20) 

This again suggest the identification of (Tyy) as the pressure P. However due to the 
breaking of the 50(3) symmetry {Txx) is not the pressure P but most also contain terms 
which are non-zero in the broken phase, i. e. terms which contain the order parameter {Ji). 
For instance it may be written as 

(r..) = P + A(Jf)(J-f), (2.21) 

where A is a measure for the breaking of the rotational symmetry and is given by 



34 [}2af-{^^fa^) ;| 

A — 2 !^ I iA (2 221 

«W u^^f (^J)2- ^'-''^ 



2 



Using this identification we can write down the stress-energy tensor for the dual field 
theory in equilibrium in a covariant form 

(r^-^) = eu^'u" + PPI''' + APf'xP%{Ja){Ja^) , (2.23) 
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Figure 1: The energy density (7tt) (a) and the charge density (1/3) (b) over the reduced temper- 
ature T/Tc for a = 0.316. The red Une is the solution without a condensate and the blue line the 
solution with (J^) ^ below T^. 

where e = (7~**) is the energy density and P^'^ = u^u^ + rj^^'^ is the projector to the space 
perpendicular to the velocity u^^. 

In figure 1 we plot {Ttt) and (^3) versus the reduced temperature, respectively. We see 
that in both cases there is one solution for temperatures above Tc and two for temperatures 
below Tc- From considerations in [22] we know that the solution with condensate (blue 
line) is the thermodynamically preferred one. For further plots see [22]. 

In addition in [22] it was found that the order of the phase transition depends on the 
ratio of the coupling constants a. For a < ac = 0.365, the phase transition is second order 
while for larger values of a the transition becomes first order. The critical temperature 
decreases as we increase the parameter a. The quantitative dependence of the critical 
temperature on the parameter a is given in figure 2. The broken phase is thermodynami- 
cally preferred in the blue and red region while in the white region the Reissner-Nordstrom 
black hole is favored. The Reissner-Nordstrom black hole is unstable in the blue region 
and the phase transition from the white to the blue region is second order. In the red 
region, the Reissner-Nordstrom black hole is still stable however the state with non-zero 
condensate is preferred. The transition from the white to the red region is first order. In 
the green region we cannot trust our numerics. At zero temperature, the data is obtained 
as described in [34,35]. 



3. Perturbations about Equilibrium 

In this section we study the response of the holographic p-wave superfluid under small 
perturbations. On the gravity side these perturbations are given by fluctuations of the 
metric hMNi-c^ jf) and the gauge field a\,j{x^ ,r). Thus we study in total 14 physical 
modes: 5 from the massless graviton in 5 dimensions and 3x3 from the massless vectors in 
five dimensions. Due to time and spatial translation invariance in the Minkowski directions. 
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0.00 0.02 0.04 0.06 0.08 

T 

Figure 2: The phase structure of the theory: In the blue and red region the broken phase is the 
thermodynamicany preferred phase while in the white region the Reissner-Nordstrom black hole is 
the ground state. In the blue region the Reissner-Nordstrom black hole is unstable and the transition 
from the white to the blue region is second order. In the red region the Reissner-Nordstrom black 
hole is still stable. The transition form the white to the red region is first order. The black dot 
determines the critical point where the order of the phase transition changes. In the green region 
we cannot trust our numerics. 

the fluctuations can be decomposed in a Fourier decomposition 

hMNix", r)= [ 4^6'^--'^ AMiv(fc^ r) , 

^ ^ J (3.1) 

To simplify notations we drop the hat on the transformed fields which we use from now on 
if not stated otherwise. 

3.1 Characterization of Fluctuations and Gauge Fixing 

In general we have to introduce two spatial momenta: one longitudinal to the condensate 
/cjl and one perpendicular to the condensate k±, i.e. A:^ = (w, fcy , A;_l, 0). Introducing the 
momentum perpendicular to condensate breaks the remaining rotational symmetry SO{2) 
down the discrete Z2 parity transformation P±: k± ^ —k± and x± — )• —x±. Thus in- 
troducing this momentum forbids the usual classification of the fluctuations in different 
helicity states of the little group since the symmetry group just consists of discrete groups 
at best P|| X P±. We do not study this case further in this paper. However a momentum 
exclusively in the direction longitudinal to the condensate or zero spatial momentum pre- 
serves the SO{2) rotational symmetry such that we can classify the fluctuations according 
to their transformation under this 50(2) symmetry (see table I). The modes of different 
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dynamical fields 


constraints 


# physical modes 


helicity 2 


hyz ) hyy hzz 


none 


2 


helicity 1 


hty^ hxy] CLy 


hyr 


4 




htz ; hxz ) Oi^ 


hzr 


4 


helicity 


httj hxxj hyy + hzz J hxtj CLt > 


htn hxr: ^rr; 


4 



Table 1: Classifications of the fluctuations according to their transformation under the little group 
SO{2). The constraints are given by the equations of motion for the fields which are set to zero due 
the fixing of the gauge freedom: a° = and hrM = 0. The number of physical modes is obtained 
by the number of dynamical fields minus the number of constraints. Due to S0{2) invariance the 
fields in the first and second line of the helicity one fields can be identified. 



helicity decouple from each other. The momentum longitudinal to the condensate, however, 
breaks the longitudinal parity invariance Py. 

In order to obtain the physical modes of the system we have to fix the gauge freedom. 
We choose a gauge where a" = and hMr = such that the equations of motion for 
these fields become constraints. These constraints fix the unphysical fluctuations in each 
helicity sector and allow only the physical modes to fluctuate. The physical modes may be 
constructed by enforcing them to be invariant under the residual gauge transformations, 
5a^ = and 5hMr = (see appendix B), 



helicity two: E = gyyhyz,hyy - h^z , 
helicity one: ^' = gyy{ijjhxy + knhty); 



(3.2) 



and helicity zero: 

^1 =^y^ 



^ ^ k^Na^f\ 2k ^ 

^3=^.-- ^2,2 ^t + -Ctx, (3.3) 
1 k , 1 , w' ■:, Sw' ^ k iuj'^w' + w(j)(j)') 



2 



u}(t)'{k'^ — w'^) 



with 



iy = gy^hyy, ^x = 9^''' hxx , 6 = 9^^ , ^tx = 9'''' ^x ■ (3.4) 



3.2 Equations of Motion, On-shell Action and Correlators 

In the following we will focus on the response exclusively due to time dependent pertur- 
bations, i.e. k^ = (w, 0,0,0). In this case in addition to the 5*0(2) symmetry, P\\ parity is 
conserved which allows us to decouple some of the physical modes in the different helicity 
blocks. In this section we write down the equations of motion for the fluctuations, deter- 
mine the on-shell action and vary the on-shell action with respect to the fluctuations to 
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obtain the retarded Green's functions G of the stress-energy tensor T^^^ and the currents 

Cf'^'P'^ik) = -iJ dM^xe-i'^^ 9{t){[T'"'{t,x),TP''{0,0)]) , 

GaW = -i / dM^xe-'^^' e{t){mt.S), J,^(0,0)]) , 

I- (3.5) 
G^'''P^{k) = -i / dtd^xe-i'^^' ^(t)([T^'^(t,x), J^(0,0)]) , 

T^'^ and J a are the full stress-energy tensor and current, respectively. Thus they include 
the equilibrium parts, {T^^) and {Ja), as well as the corresponding dissipative parts which 
arise due to the inclusion of fluctuations in our model. In the following we split the analysis 
into the different helicity blocks. 

3.2.1 Helicity two mode 

First we look at the non-trivial helicity two mode displayed in table 1. If we expand the 
action (2.1) up to second order in the fluctuations, this mode decouples from every other 
field. Therefore it can be written as a minimal coupled scalar with the equation of motion 

The contribution from this mode to the on-shell action is 




'helicity 2 ^2 J (2vr 



i_ r 

2 
5 



~ 2f~ W~ 2^)^ ~ 



=^bdy 

(3.7) 

which is divergent as we send rbdy oo. The divergence can be cured by holographic 
renormalization (see appendix A). The renormalized on-shell action is 

.4 , ^4^ 



oon-shell 

•^helicity 2 - ^2 y (27r)4 



-0^2 



2\ " "^32 



(3.8) 



where cu = w/r/^ is the dimensionless frequency, Hq and are defined similarly to the 
quantities in (2.11) and SqS2 is a short form for Hq(w)S2(~'^)- Now we use the recipe by 
Son and Starinets [36] to compute the Green's function of this component. The response 
due to the perturbation hyz is given by 

^^on-shell 

with 

-(fifM-<^"> + ^'^0' <^-'"' 

where {Tyy) is the equilibrium contribution given by the pressure P. As we will see in section 
4.1, the Green's function of this helicity mode will lead to a shear viscosity component with 
universal behavior, i.e. r/y^/s = l/47r. 
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3.2.2 Helicity one modes 



Now we look at the helicity one modes displayed in table 1. Again we obtain their equations 
of motion by expanding the action (2.1) up to second order in the fluctuations and varying 
it with respect to the corresponding fields. The equations of motion are 

^ 3// (l 2f' N' a'\ 3/ / oj^ /W 2a^cP'^\ 3 ^ , 



and 



"^t + ^al (3.11b) 



+ „2^27\r2_2% „2^2Ar2_2%' {6.lZa.) 



1" , 2/' N' a'\ 1/ g , , f UJ ^ ^^ I 



,2 



'^"'^.a^, (3.12b) 



=a,, + — + — ^ a„ + „ + „ a. 



where = g^^hty and ^'x- = g^^hxy Note that due to the parity Py, the helicity one 
modes split into two blocks where the modes of the first block are even while the modes 
of the second block are odd under Py. In the first block there is only one physical mode 
Uy while the value of the other field ^'4 is given by the constraint (3.11b). This can also 
be seen in the gauge invariant fields (3.2) since hty drop out for k^i^ = 0. The other three 
physical modes appear in the second block where = ^ for k^^ = 0. 

The contribution from these modes to the on-shell action is 



^on-shcU 
'helicity 1 



1 /-d^fc fr^^ 1 
(27r)4 \ 4ct 4 



2a 



+ 



1 



2 r-^pNa 



N 

+ —4 — 4*- 



2a 



at^t 



N 



=^bdy 



ra^Na 

2 2f' N' 
~ 2N 



11/ 22' 3*^' 



a 



a 



(3.13) 
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which is again divergent^. The renormahzed on-shell action is given by 



-lon-shell 
'helicity 1 



b 



/-2\^/-~2\^ 

K)oK)i + K)oK)i + K)oK)i 



62 



6 ~f, /~1\f)/,T. \b' 



(3.14) 



where = r/^aj^ is dimensionlcss. We obtain the response of the system due to the 
fluctuations and hfy by variation of the on-shell action, 



I 



f c-on-shcll \ 
"'^helicity 1 \ 

r con-shcll 
"'-^helicity 1 



'Gl 



3,3 



,0; 




(3.15) 



with 



^G|;|(a;) Gr^(a;) 
,G*f|(a;) G*f'*f(a;) 











HI 


'7 ( 







2(ag);H 




(3.16) 



This result agrees with the result obtain in the holographic s-wave superfluids [17,30] and 
thus the breaking of the rotational symmetry has no effect on this subset of fluctuations. 
The coupling between the current and the momentum Tty is known as the thermoelectric 
effect which we will study in the next section. 

The response due to the fluctuations a^, and hxy is given by 



hclicity 1 

^Con-shell 

helicity 1 

5(^r 



Si con-shell 
^'^helicity 1 



lG\l{^) Gl^iu:) Gr{u)\ 

Gfyiu) Gfy{u) Gr{u) 
vG^f^(a;) G^yl{u) G'^V'^^y {oj) j 



where the matrix of the Green's functions is given by 



w 

h 

'^5 



2 Mm 

6 



'^5 



3 



2^^) , 1 .-,4 
+ 32^^ 



(3.17) 



(3.18) 



^Note that the contribution of the on-shell action is zero at the horizon since we can set *t to zero there. 
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Due to the breaking of the rotational symmetry we see a new couphng between the currents 
J\2 a-iid the stress tensor T^y in this subset of the fluctuations. This new coupling generates 
some interesting new physical effect: it induces a non- universal behavior of the ratio of shear 
viscosity to entropy density and a flexoelectric effect known from nematic crystals. We will 
study these effects in the next section. 

4. Transport Properties 

In this section we extract the transport properties of the holographic p-wave superfluid 
from the correlation functions presented in the previous section. We split our analysis into 
distinct transport phenomena. 

4.1 Universal Shear Viscosity 

Let us start by considering the helicity two mode hy^. It is well known that, in the isotropic 
case, the corresponding component of the energy-momentum tensor may be written as^ 

{Ty') = -{P + \ujr]y,)hy,. (4.1) 

Using (D.4) we see that this result is still correct in the transversely symmetric case we 
are studying here. The result also agrees with our gravity calculation (3.10) ^ . Thus the 
shear viscosity is given by the well-known Kubo formula 

^ = - hm -Im {Cy^'y^) = - hm A ^ ^4 2) 

In the following we show that we can apply the proof for the universal result for the ratio 
of the shear viscosity to entropy density described in [25]. In the uj ^ {) limit the equation 
of motion (3.6) corresponds to (9^11 = 0, with 11 the conjugate momentum to the field 
S = gyyhyz- This is the decisive condition in order to apply the proof of [25]. Therefore we 
conclude that here we obtain the universal result for the ratio of shear viscosity to entropy 
density [23-25], 

— = • (4.3) 
s All 

In this subset we do not see any effect of the rotational symmetry breaking since the 
fluctuation hyz is transverse to the condensate. 

4.2 Thermoelectric Effect perpendicular to the Condensate 

Now we relate the results of (3.15) to the thermoelectric effect on the field theory side. We 
begin with the well known connection between electric (J3") = {J'^) and thermal (Q"*") = 
(T*-*-) — /u(J-'-) transport perpendicular to the condensate direction (see e.g. [17,30,37] for 
the same calculation for holographic s-wave superfluids), i.e. 

(4.4) 

^Note that g"" = j?^" - /i^". 

^We do not see an a;*-term as in (3.10) in the Unear hydrodynamic description since this term corresponds 
to higher order term with four derivatives. 
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where the electric field E± and the thermal gradient —'V±T/T need to be related to the 
background values of the gauge field (aj_)Q and the metric (^t)o- This identification is 
done in [30], 



(4.5) 



Putting all together and comparing the relation of the electric and thermal transport to 
the corresponding equations for (J^) and (T*^) in (3.15), we can identify the transport 
matrix of (4.4), 



±± 



a 



Ta 



±± 



TR 



±± 




±± 



(4.6) 



{Ttt) + fi^a 



2_±± 



These results are in agreement with [30]. The coupling between thermal and electrical 
transport is well known in condensed matter physics, since the charge carriers (electron or 
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Figure 3: Real part of the conductivity Re((T^^) over the frequency lo/{2ttT) for a = 0.032. The 
color coding is as follows: blue T = oo, red T = 1.34Tc, brown T = LOOT^, green T = 0.40T;, 
orange T = O.lOTc- Note that the three curves with the highest temperature, blue, red and brown, 
are nearly on top of each other. The agreement of the curves in the cj — > limit is due to the small 
change in the strength of the Drude peak with temperature. Below Tc, the supcrfluid contribution 
to the delta peak at w = is turned on and we obtain larger deviations from the T ~ oo curve, 
since the area below the curves has to be the same for all temperatures (sum rule). Furthermore 
the value for w — >■ clearly asymptotes to with decreasing temperature. 
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Figure 4: Real part of the conductivity Re((T^^) over the frequency lli/{2ttT) for a = 0.316. The 
color coding is as follows: blue T = oo, red T = l.OOTc, brown T O.SSTc, green T = 0.50Tc, 
orange T = 0.19Tc. In this plot we see that the Drude peak has already a much stronger dependence 
on the temperature than in the a = 0.032 case, since the blue and the red curve can be clearly 
distinguished. Below Tc the contributions of the superfluid phase to the delta peak leads again to 
a tendency of the curve to vanish for frequencies in the gap region, since the area below the curves 
have to be the same for all temperatures (sum rule). 

holes) transport charge as well as heat. In this subset we do not observe any effect of the 
breaking of the rotational symmetry since all the fields are in the transverse direction to 
the condensate. 

In figure 3, 4 and 5 we plot our numerical results for Re(<T"'""'") versus the frequency 
uj/{2ttT) for different values of a as defined in (2.3), namely a = 0.032 < Oc, a = 0.316 < etc 
and a = 0.447 > Uc-, respectively. For large frequencies, i. e. w ^ 27rT, the conductivity 
asymptotically has a linear dependence on the frequency (e.g. [38]), 

I I 

Re(a-^-^) ^ ^TTCj for u 2tiT . (4.7) 

For small temperatures, i. e. T < Tc, we see a gap opening up at small frequencies. The 
size of the gap increases as the temperature is decreased. This is the expected energy gap 
we know from superconductors. The gap ends at a frequency ujg with a sharp increase 
of the conductivity. Beyond the gap the conductivity at small temperature, i. e. T < Tc, 
is larger than the corresponding value at large temperature, i. e. T > Tc such that the 
small temperature conductivities approach the asymptotic behavior (4.7) from above (cf. 
[18,19,21]). 

The value of Re(<T"'-"'") at a; = approaches zero with decreasing temperature. Below 
Tc the tendency for this part of the conductivity to vanish increases. Nevertheless, we still 
find finite values even below Tc, i.e. these values seem to be suppressed but not identically 
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Figure 5: Real part of the conductivity Re((T^^) over the frequency lli/{2ttT) for a = 0.447. The 
color coding is as follows: blue T = oo, brown T = 1.95Tc, green T l.OOTc, red T = 0.91T;, 
orange T = 0.34Tc. Again we see the same tendency as before for the curve to vanish at a; ^ 
for decreasing temperatures. The strength of the Drude peak has a strong dependence on the 
temperatures, since the blue and the brown curve are quite far apart (both curves were computed 
for temperatures above T^). 

vanishing (c.f. [17]). In [34] it is shown that in the limit T — ?• there is a hard gap, i. e. the 
value for the conductivity becomes zero. Finally, we observe that an increase in a leads to 
a stronger suppression of the real part of the conductivity in the gap region. 

Due to the sum rule for the conductivity, i.e. the frequency integral over the real 
part of the conductivity is constant for all temperatures, a delta peak has to form at 
zero frequency which contains the "missing area" of the gap region. The strength of the 
delta peak has two contributions: the first is proportional to the superfluid density n^, 
Re(o"-'-'-) ~ oP' / 1^1 TTns6{uj) and appears only for temperatures below Tc. The second 
contribution is a consequence of translation invariance of our system, the Drude peak, and 
appears for all temperatures. 

The delta peak is observed in the imaginary part of the conductivity by using the 
Kramers- Kronig relation (see [17]), 



Drude peak. In figure 6 we present the imaginary part of the conductivity a;Im(fT ) versus 
the frequency oj/{2t^T) for a = 0.316 and different temperatures. We see that (jjlm.{a^^) 
takes finite values for a; — )• and T < oo. The finite values above Tc are due to the Drude 
peak, i.e. the part of (4.8). Below we see a further contribution from the superfluid 




(4.8) 




parametrizing the contribution from the 
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Figure 6: Imaginary part of the conductivity a;Im(cr^''-) over the frequency uj/{2'kT) for a — 0.316. 
The color coding is as fohows: blue T = ex., red T = LOOT^, brown T = 0.88Tc, green T = 0.50Tc, 
orange T = O.lOTc- The curves in this plot have a constant value for 2^ ~^ 0, which is determined 
by the (S-peak in the real part of the conductivity by the Kramers-Kronig relation. The values for 
this constant are, in the same order as the temperatures above: 0, 8.0, 14.7, 50.3 and 362.4. Note 
that we already see a finite value for T = l.OOTc, this is due to the Drude peak. Below Tc the 
values at w = are due to two contributions, first the Drude peak, as before, and second due to 
the superfluid density. 

density. By analyzing the temperature dependence of Ywii^^Q uilvii{a^^) , we get a smooth 
curve, which is, however, not differentiable at for a < Oc, i.e. it behaves as equation 
(4.8) anticipated. However, for a > etc we see a jump at Tc as consequence of the jump in 
the condensate. Furthermore Ad and Ag have a non trivial dependence on a. Finally, as 
expected, there is an increase in the superfluid density with decreasing temperature. 

4.3 Non-Universal Shear Viscosity and Flexoelectric Effect 

Now let us study the remaining three components of the helicity one modes, (Jj*"), {J2) 
and (T^-*-) as given by (3.17). We first focus on (T^"*"), which for (a^)^ , (0^)0 = can be 
translated into the following dual field theory behavior 



where r]x± is the second shear viscosity which is present in a transversal isotropic fluid (see 
appendix D) and with {Txx) as defined in (2.19). Here we see again that we can apply the 
Kubo formula to determine the shear viscosity r]x±, 



As described in [28], this shear viscosity has a non-trivial temperature dependence even 
in the large N and large A limit and is therefore not universal. In fig. 7 we compare our 
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(4.10) 
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Figure 7: Ratio of shear viscosities rjyz and ?7^j^ to entropy density s over the reduced temperature 
T/Tc for different values of the ratio of the gravitational coupling constant to the Yang-Mills coupling 
constant a. The color coding is as follows: In yellow, rjyz/s for all values of a; while the curve 
for rix±/s is plotted in green for a — 0.032, red for a = 0.224 and blue for a = 0.316. The shear 
viscosities coincide and are universal in the normal phase T > Tc- However in the superfluid phase 
T < Tc, the shear viscosity r]yz has the usual universal behavior while the shear viscosity rix± is 
non-universal. 

numerical results for the ratio of the shear viscosity r]x± to the entropy density s with the 
universal behavior of the shear viscosity r]yz for different values of a. We see that in the 
normal phase T > Tc, the two shear viscosities coincide as required in an isotropic fluid. 
In addition, the ratio of shear viscosity to entropy density is universal. In the superfluid 
phase T < Tc, the two shear viscosities deviate from each other and rix± is non- universal. 
However it is exciting that r]x±/s > l/47r, such that the KSS bound on the ratio of shear 
viscosity to entropy density [23] is still valid. 

The difference between the two viscosities in the superfluid phase is controlled by a 
as defined in (2.3). In the probe limit where q = 0, the shear viscosities also coincide in 
the superfluid phase. By increasing the back-reaction of the gauge fields, i. e. rising a, the 
deviation between the shear viscosities becomes larger in the superfluid phase as shown 
in fig. 7. If a is larger than the critical value ac = 0.365 found in [22] (see fig. 2) where 
the phase transition to the superfluid phase becomes first order, the shear viscosities are 
also multivalued close to the phase transition as seen in fig. 8. Since there is a maximal a 
denoted by amax = 0.628 for which the superfluid phase exists (see fig. 2), we expect that 
the deviation of the shear viscosity rix± from its universal value is maximal for this amax. 
Unfortunately numerical calculations for large values of a are very challenging such that 
we cannot present satisfying numerical data for this region. It is interesting that also the 
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Figure 8: Ratio of sliear viscosities rjyz (blue) and r]x± (red) to entropy density s over the reduced 
temperature T/T^ for a = 0.447, which is larger than the critical value where the phase transition 
becomes first order: The shear viscosities coincide in the normal phase T > and are universal. 
In the superfluid phase r]xi^ is non-universal. Close to the phase transition, it is multivalued as 
expected for a first order phase transition. 

deviations due to A and Nc corrections are bounded. In this case the bound is determined 
by causality [39]. 

As described in our letter [28], we also have found numerically that for a < Oc 



Interestingly, the value of /3 appears to be independent of a. This result has recently been 
confirmed by an analytic calculation in [40]. 

The non-universality of the shear viscosity can be understood in the following way. 
For the r]x± component, the equation of motion (3.12a) in the w — t- limit includes also 
non- vanishing source terms besides the derivative of the conjugate momentum Hx of 
i.e. dr^x = source. This is in contrast to equation of motion which leads to the r]yz 
component. Note that the source term depends on the condensate w and the fluctuation 
a\ and vanishes if the condensate w vanishes. Hence, as we confirm numerically in fig. 7, 
when the condensate is absent (i.e. for the T > Tc case) we obtain again the universal 
result, since in this case the same proof as described above for the helicity two mode applies. 

1 2 

For hx± = 0, we obtain flavor charge transport, i. e. a flavor field generates a 
flavor current {J 12)- In the unbroken phase it is useful to combine the fields a^^ in the 
way = lb \a\ since they transform in the fundamental representation of the U{l)z 
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symmetry and are complex conjugate to each other. To make contact to the unbroken 
phase, we also use this definition in the broken phase. 

We also use the definition for the currents (J±) = 1/2 ((Ji) ± i( J2)), such that the full 
transport matrix becomes 







/ 
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\ 



QX±- 
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\hx±/ 



(4.12) 



where the flavor conductivities are given by 



G 



Gi;^(a;) 



QXA 



± 



.CO 



G: 



x± , 




Gi^^{u)±iGi 



(4.13) 



First note that for /x = where the SU{2) symmetry is restored, i. e. = a^, the Green's 
function is diagonal and G^''\i = G't''^ = G^'^. In the unbroken phase, G'\^''t = G'^''^ = 
is still valid for /x 7^ 0, since the do not couple to each other, while Glji'l}: 7^ Gl'i'i'l for 
// 7^ 0. In the unbroken as well as in the broken phase we find 



G^l^uj) = G'^'^^i-uT , Gi;^(a;) = G^;i:(-^)* 

Gi^'^iuj) = G^^^{-ujy and G''^+{u) = G^^U-^) 



(4.14) 



as expected since a^{u}) = (a"^(— w))*, a^(cj) = (a^(— w))* and ^x(w) = {^x{—co))*. 

In figure 9 we plot the real and imaginary parts of G_|_'_|_((x!). We see in fig. 9(c) and 9(d), 
showing Im(G^'^), that for temperatures T > Tc, the quasinormal modes tend towards the 
origin (in these plots we see their projection on the real axis) (see e. g. [18,41]). For T <Tc 
we see a pole at the origin which is due to the massless Nambu-Goldstone modes. These 
Nambu-Goldstone modes are related to rotations of the director (J^^) in real space which 
are generated by the fluctuations a\ ^. Furthermore, as expected for large frequencies. 



*The other Nambu-Goldstone mode is related to the change of the phase of the condensate and corre- 
spond to the fluctuation which shows up in the helicity zero sector. 
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Figure 9: These plots show the real and imaginary part of the correlators G^'^ versus the reduced 
frequency ijj/{2ttT) for a = 0.316 at different temperatures: T = oo blue line, T = 3.02rc red line, 
T = l.OOTc brown line, T = OMT^ green line and T = 0.50Tc orange Hue. 




(a) (b) 

Figure 10: These plots show the real and imaginary part of the correlator G'^''^ versus the reduced 
frequency uj/{2ttT) for a — 0.316 at different temperatures: T = oo blue line, T — 3.02rc red line, 
T = l.OOTc brown Hue, T = O.SST^ green hue and T = O.SOTc orange line. The curves for the 
temperatures above Tc are exactly zero for all frequencies. 



the Green's function grows proportional to in the zbib components as for the correlator 
Ggg"^. In figure 9 the correlators for the different temperatures do not seem to have the 
same asymptotic behavior. However, in the present case we have contribution from terms 
such as w/i, i.e. of first order in a;, which are not existent in the G^'^ component. Hence to 
see that all correlators have the same limit, larger values of uj have to be considered. Even 
not present in our figures we verified numerically that the asymptotics of the correlators at 
different temperatures agree. Furthermore, in this context, we check that the fluctuations 
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Figure 11: These plots show the real and imaginary part of the correlators G±. 
reduced frequency uj/{2ttT) for a = 0.316 at different temperatures: T = oo blue line, T = 3.02rc 
red line, T = l.OOTc brown line, T = 0.88Tc green hne and T = O.SOTc orange line. The curves for 
the temperatures above Tc are exactly zero for all frequencies. 



aj_, which are unstable in the normal phase, are stabilized in the broken phase for T < 
i.e. the quasinormal modes of these perturbations stay in the lower half plane. Thus, the 
preferred direction induced by the current (v7f ) is stable^. A more detailed study of this 
sector in the probe approximation is in preparation [42]. 

In figure 10 we plot G'ljl'l'l(ci;). We see that this correlator vanishes since the do 
not couple in the unbroken phase. Furthermore, below Tc, a pole at a; = due to the 
Nambu-Goldstone mode appears. We do not show Gl'l'ljl(a;) since Glji'i'l and G'^l'j^ look 
alike. Nevertheless, there is a difference between them in the broken phase. The difference 
arises from the contributions to the correlators due to the mixed terms, G^'^ and G^^"^, in 
the corresponding equation in (4.13). However, these are suppressed in relation to G^'^, 
which contains the Goldstone mode in the broken phase. 

Finally, in figure 11 we show G^ (uj). Note that the contribution of (J^f) to the real 
parts of the correlators is not included in the corresponding plots since it just shifts the 
curves by a constant. Furthermore, we have that Im(G^^ ) = Im(G^"'"_|_) and Re(G^^ ) — 
const = Re(G^"'"^), i.e. there is a constant offset between the real parts of these correlators. 
We expect that this constant offset is generated by the term (^xY2^^^/ ('^y)o('^) '^hich may 
be constant in the limit a; — )■ since has a normal mode at w = and the subleading 

^Note that the Coleman-Mermin-Wagner-Hohenberg theorem does not apply here since the field theory 
is in 3+1 dimensions. 
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term of is probably not sourced. A analytic calculation similar to [40] may confirm this 
claim. In the unbroken phase these correlators vanish since the differential equations of the 
corresponding fields decouple. In the broken phase the correlators present a rich structure, 
which we cannot fully address at present. However, it seems that the coupling between the 
flavor fields and the strain hx± generates new quasiparticles which appear as bumps in 
the curves. 

In addition to the flavor conductivity and the shear viscosity we obtain a coupling 
between the stress (T^-*-) and the flavor fields as well as the currents ( J^) and the strain 
hx± described in (4.12). This coupling introduces an effect which is called flexoelectric effect 
in nematic crystals [14] and only appears in fluids with broken rotational symmetry. We 
have a current {Ji) in a favored direction in the background which interacts with the 
flavor flelds a^. This interaction induces a force on the current which pushes the current 
in its perpendicular direction generating the stress (T^"*"). In the similar way, a strain 
hx± introduces an inhomogeneity in the current ) resulting in a flavor field which 
generates the currents {J±). 

5. Conclusion 

In this paper we have studied transport phenomena in holographic p-wave superfluids con- 
structed in the SU{2) Einstein- Yang-Mills theory. We classify the perturbations about 
equilibrium according to their transformation properties under the symmetry group. At 
zero momentum, there is an SO{2) symmetry left which allows us to divide the perturba- 
tions into different helicity sectors: helicity two, one and zero states. While the helicity 
two state is trivial and leads to the universal ratio of shear viscosity to entropy density, the 
helicity one states are non-trivial. Due to a Z2 parity, this sector splits into two blocks. In 
the first block we find the thermoelectric effect transversal to the direction favored by the 
condensate. In the second block we obtain two interesting new phenomena: a non-universal 
shear viscosity and a flexoelectric effect. These two effects are due to the anisotropy of our 
system. 

Anisotropic fluids have been studied in particular in the context of nematic crystals 
whose hydrodynamic description is given in [14,43]. In this paper, we have initiated 
the connection of the hydrodynamic description of anisotropic fluids with gauge/gravity 
duality. The results we obtain in this paper are in agreement with this description, i. e. the 
transport coefficients found here can be related to the ones in [43]. However since we have 
not studied the helicity zero modes in much detail, we have not yet described all transport 
properties. In particular, the thermoelectric effect along the condensate as well as the 
coefficients (x, Cy and A described in appendix D are still missing. In the future we plan to 
study these coefficients in detail. This study may also lead to a covariant hydrodynamic 
description of anisotropic superfluids. 

Furthermore, analytic results close to the phase transition can be found for small values 
of a [40], which determines the ratio of the gravitational to Yang-Mills coupling. On the 
one hand this analytic approach allows for a detailed study of the transport coefficients 
close to the phase transition. On the other hand it also permits us to use the fiuid/gravity 
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correspondence in order to obtain the complete hydrodynamic description of the system 
directly from gravity. Similar analyzes for holographic s-wave superfluids can be found 
in [10-12]. We intend to follow this line of thought further. 
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A. Holographic Renormalization 

The goal we are pursuing in this section is to find covariant counterterms which can be 
subtracted from the action (3.7) and (3.13) in order to make it finite. We follow the lead 
of the references [44, 45] to perform the holographic renormalization. 

A.l Asymptotic Behavior 

In this section we look at the behavior of the fields at the horizon and at the boundary. We 
want to calculate real-time retarded Green's functions [36,46], i.e. at the horizon, besides 
regularity^, we have to fulfill the incoming boundary condition. For this purpose we plug 
in the ansatz, 

i^(r)U^ = 6fE4i^'' ^ith e. = -^-l, (A.l) 

for the behavior of the fields near the horizon, into the equations of motion (3.6), (3.11) 
and (3.12). It turns out that, as expected, we obtain two possibilities for /3, namely 

/3 = ±«7^, (A.2) 

with T being the temperature defined in equation (2.12). As said before, we choose the 
solution with the "— " sign which corresponds to the incoming boundary condition. Note 
that the other solution represents the outgoing boundary condition. 

Our ansatz at the boundary is similar to the one used for the background calculation 
in section 2. However, here we have to add a logarithmic term to get a consistent solution 
(c.f. [44]). Therefore we use 

^(^)U.. =E4(^' + ^^'lne,) , with 6,= (^)'. (A.3) 



j>0 



r 



Let us now use the above expansions for the helicity one states. In the case of the 
equations (3.11) we have 3 independent expansion coefficients at the boundary (4 free 
parameters from the 2 second order differential equations minus 1 free parameter due to 



®Even with all fluctuations switched on, there is no need for a further constraint besides cpirn) = at 
the horizon to guarantee regularity. 
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the constraint). We choose them to be (a^)^ , [o-y)^ and (^t)Q- At the horizon we aheady 
halved the independent parameters by choosing the incoming boundary condition. From 
the remaining two parameters we can get rid of one by using the constraint equation 
(3.11b). Therefore we are left with one free parameter at the horizon, we choose (ay)g- 

When solving these equations numerically we set {ciy)^ = 1 and scan through different 
values of uj. 

We can perform similar considerations for the second part (equations (3.12)). In this 
case it is even simpler. We do not have any constraint, just three fields and their corre- 
sponding equations of motion. Therefore at the boundary we have six independent pa- 
rameters, namely {al)^, i^Do^ (^^Oo ^^"^ ('^^■)2' horizon we have 
(^a;)g, (a^)g , {o-Dq- Note that, as before, we already fixed three free parameters at the 
horizon by choosing the incoming boundary condition. Again by choosing the values for 
all fields at the horizon the system is fully determined. 

Notice that the same is true for the helicity 2 state. We have again 2 independent 
components at the boundary, namely (S)g and (S)^ which are fixed by the incoming 
boundary condition and (")g at the horizon. Therefore as before, the equation is fully 
determined. 

In the following we state the first few non-vanishing coefficients of the expansion at 
the boundary of the different fields. We need them later on to determine divergences in the 
on-shell action and to calculate the Green's function. The explicit form of the dependent 
coefficients are 



1 r 
2 

1 r 
2 

i4 



UJ {ay}^-2iuj^i{ay}^ + fi (aj^ 
h ' 

'''' ' -4/^n6 



-a 



21b ^~3\f> 



16 

b 



3 

b _ 1 
2 ~ 16' 



(A.4a) 
(A.4b) 
(A.4c) 
(A.4d) 
(A.4e) 
(A.4f) 



Note that fi = (pQ and (pi are the expansion coefficients of (/){r) at the boundary. 

We do not state the expansion at the horizon, since there is no additional information 
to equation (A.l), and the explicit form of the non-independent coefficients is very long. 
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A. 2 Counterterms 

By plugging the expansions (A. 3) into (3.7) and (3.13), resulting in 



mr 



m, 



fe2 



(4)o(4)t + («^)o(«y)t + («y)o(«y)l 



1 



2-2 

a /i 



~b/~l\b 



^(aJ);(^.)^ + 2c^^(a3)J(^, 



+ \a 



''^A(«y)o(«y) 



(1 -2 "^-4 "^~4i 

+ ( offo'^ "^64^ "^32^ 



1 2-2 2~2i 

-a u — -a u hi eh 



i^lYo + i^lYo' + i^X 



(A.5) 



2 ~ 2 1 

+ -a Inefo 



a: 



^'^A(«y)o(«y)olne6 



la 



=^bdy 



Note that Son-shcii = J '^:Jyi ^ri, with r = rbdy ^ 1- Moreover, as before, the first field 

expansion coefficient is always a field of —k and the second of k, e.g. (H)p(— /c) (H)2(fc). It 
is obvious from the last few lines of (A.5) that we have to add counter terms to the on-shell 
action S'on-shcii to take the divergences in r into account. 

The terms that have to be considered are the ones in (A.5) with explicit r dependence. 



1 -2^ 1 ~4i 



Ho +(^-)o 



1 



-a In Eh {Co 



,-.2 



Yo + i^lYo + i^lYo 



and 



la 



First we need the induced metric 7 on the r = rbdy plane. The induced metric is defined 

by 



dx 



(A.6) 



=^bdy 



resulting in 



-A(rbdy)cT(rbdy)^dt2 + 



Thdy 



bdy 



/(^^bdy)' 



^dx2 + r^dy/(^bdy)^(dy2 + dz^). 



(A.7) 



We do not literally derive the covariant counter terms here in this work. However, by 
looking at the counter terms of B. Sahoo and H.-U. Yee calculated in [45] we get an idea 
how they should look like, namely some combinations of -^[7], R^vY^] and F^^. The first 
two are the Ricci scalar and Ricci tensor on the induced surface respectively, the latter 
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is the field strength tensor on that surface. Possible covariant combinations of the three 
terms are ^J—'^R[^], ^—^W^ [7] R^j^y [7] and ^ —^F^^F'^^^ . Now lets have a look at their 
expansion for r ^ 1. 



^R^-Yi\R^,M 



'0 



fe2 



and 



2 



K)o' + («^)o' + (4)f 



vb2 



«Do +(4)o -4i^/^(4)o(4) 







Therefore by adding the real space action 



1 



d^x 



\Rb] - ^R^^^b]R,M In £6 + yF;,F"^^ln6,) 



(A.9) 



'bdy 



to the action Son-shell (3-7) and (3.13) we get a divergence- free theory (up to the second order 
in the fluctuations) for rbdy S> 1, i.e. also the real time Green's functions are divergence 
free. The renormalized rbdy ^ 1 Lagrangian is then given in (3.8) and (3.14). 

B. Constructing the Gauge Invariant Fields 
B.l Residual Gauge Transformations 

The transformations we have to look at are diffeomorphisms and SU{2) gauge transforma- 
tions. On the one hand, we demand that the fields be diffeomorphism invariant, i.e. 

Sj^^ = Ct.^ = 0. 

£s is the Lie derivative along S, i.e. 

C-T,gMN = ^ M^N + N^M = Om^N + Qn^M — '^^AIN^P' 

with r£fjy being the Christoffel symbol. 

On the other hand they have to be invariant under the SU{2), 

$ ^ M(A)$ = 

with M(A) being the SU{2) transformation matrices, this is equivalent to 

5a ^> = 0. 



(B.l) 



(B.2) 



(B.3) 



(B.4) 

$ stands for the physical modes in our system and is composed of the helicity fields 
^txi £,y, Ox ^nd of, with a = 1,2, 3. The invariance of $ under the above transforma- 

tions translates into 

3 

<5$ = (5s + 5a)$ = 5^(raK + ) + Tjd^tx + TgS^t + TgS^x + mS^y = 0, (B.5) 

a=l 

with Ti being the r dependent coefficients. 
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r** = 

rr 


C4' 


^ tr 


pr- _ 


ClCl' 





B.1.1 Diffeomorphism Invariance 

Let us look at the invariance under diffeomorphisms. We begin by defining 

guN = 9MN + hMN, ,^ ^, 

(B.6) 

— + (^M- 

Furthermore note that Sm (and later on A) are of the same order as the fluctuations. 
Through the gauge choice hMr = we can determine the form of Sm up to some constants, 
because 

c^mr = ^ ^^^^ 

=>^M5^r + drT>M — 2^J^^^Sp = 0. 

Note that we just need the Christoffel symbols to zeroth order in fluctuations, i.e. the 
background Christoffel symbols. They are 

' xr ^ vr ' 

, , (B.8) 

C2C2 , rr - '^3*^3 

2 ^ VV ~ 2 ' 

C4^ C4^ 

with 

ds^ = QMNdx^dy^ = -ciirfdf + C2{rfdx^ + C3{rf{dy^ + dz^) + c^irfdr^. (B.9) 

We get 4 equations (+1 for the z component which is exactly the same as the one for the 
y component), which read 

- iij^r + St' - 2— St = 0, (B.lOa) 
ci 

ikEr + J^x -2—Y.^ = 0, (B.lOb) 

Sy'-2^Sj/ = 0, (B.lOc) 

2E/-2— Er = 0. (B.lOd) 

C4 

We work in momentum space, i.e. the ansatz used is 

EM(t,x,r) = y"dS e*^'^^''S(w,fe,r), (B.ll) 
with k^^ = {oj, k, 0, 0). The solutions to the equations above are 

T,ti<^,k,r) = Ktcf + iuKrCi^A, with A = J dr- 

Ylx{uJ,k,r) = KxC2 — ikKrC2^B, with B = J dr-^, ^-g j^2) 

1:^(0;, k,r) = KyC3^, 

T,r(u!, k, r) = KrCi, 



C4_ 
2 ■ 



Cl 
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with Ki being constants. Using these solutions in the remaining equations fe^j, we get 





2iu 






















= 0, 






— 0-^x 




C2^ 



(B.13) 

Here we see already that is a physical mode! 

Applying the same procedure to the gauge fields, i.e. 5sa)i = ^E^^, it results in 

OS«x — — 2 2"^' 

3 <?> ^ , ^^^0^ 
OEflj = — ^^r H ^^t- 

The Lie derivatives of the remaining components vanish. 



B.1.2 SU{2) Gauge Invariance 

This transformation only affects the gauge fields, therefore we do not have to care about 
the metric fluctuations here. A field in the adjoint SU (2) representation transform under 
the SU{2) as 

5Ka% = VMA^it, X, r) + e'^^'^^^.A^ (B.15) 
Again we constrain the possible A" by using the gauge choice = 0, i.e. 



= 6aK = V,.A"(t,x,r) = drA''{t,x,r) 
^A''{t,x,r) = A%t,x). 



(B.16) 



We choose the ansatz A°-{t,x) = J d^x e*^'^^^A"(a;, A;) with A:^ = (a;,A;,0, 0) to calculate 
6\a'^. Note that by using the definition above of A;^ the y and z components do not mix 
with the rest and we can forget about them. We end up with 

d^al. = i/cA\ 6\a\ = -\loA^ - cpA"^, 

SAal = ikA^ - wA^, dxal = -iujA^ + cfiA^ , (B.17) 
d^al. = ikA^ + wA^, 6Aaf = -iuA^. 

B.2 Physical Fields 

Plugging everything into equation (B.5) results in 6 equations, due to the fact that S^, S^,, S^, A^, A^ 
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and are independent. The equations are 



=ikTi — iu}T4 + 0r5, 

=ikT2 + WTs — 4)Ti — ia;r5 , 

= — WT2 + \kT^ — iwTg , 

= - + 2ikTg + ikwTi, (B.18) 
=2iu}Ts + i/c^ry + iuj(j)Te, 
=w'ti + 4>'t^. 

The four physical fields, $i = Tj-(Helicity fields), we get by solving above equations, 
are 



$1 

(p {k^-w^)4> {k^-w^)4> 

^ ^ k^ci\ 2k. 



1 /c 1 1 ^ w' ct>w' ^ k iuj'^w' + W(j)(f)') 
$4 =a^ + -a^ - -w^x - -jat + — 77^2 2w7" 



3 



C. Numerical Evaluation of Green's Functions 



Here we review the algorithm to evaluate Green's function when there is operator mixing 
[47]. For concreteness we present the algorithm for the second block of the helicity one 
fields where we have the mixing of the fluctuations 0^,0^,^' 3.. Let us first go back to the 
action (3.13) and have a look at the components we need. 



2/2 



ra^f^Naw' , ^ 
+ —h, al^x 



= »-bdy 



d^k 

(2^ 



[A^{-k, r)A{k, r)drh{k, r) + A^i-k, r)B{k, r)A{k, r)) 



='"bdy 

(C.l) 



A'^{k, r) = [ay{k, r),ay{k, r), "^x{k, t)) , A{k, r) and B{k, r) are the matrices containing the 
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coefficients of the field's bilinear forms. Their explicit form for r = rbdy ^ 1 is 



/ -a^w^rl \ 
. (C.2) 

Note that all components of B{k,r\^^y) are contact terms. 

The next step is to define the boundary condition at the horizon for the vector A{k, rh). 
This includes the incoming boundary condition we already introduced before (c.f. eq. 
(A.l)). We define 

Ha)ik, t^th)^ (ef,) + O (e,)) , (C.3) 

the index / refers to the three fields a^, and '^x- Due to the fact that we have 3 coupled 
second order differential equations we need 6 independent boundary conditions. We already 
halved them by demanding incoming boundary conditions. We determine the remaining 3 
conditions by choosing 3 linear independent vectors e(i), 6(2) and 6(3), with 

6^) = (1,0,0), 6^) = (0,1,0), e5) = (0,0,1). (C.4) 

Note that alternate choices are possible, but we get the best numerical result (with least 
noise) using these values. We are now able to generate three linearly independent solutions 
for the equations of motion by solving them numerically using the three linearly inde- 
pendent boundary conditions. We put them together in a 3 by 3 matrix H{k, r) defined 

by 

H',{k,z) = Kl~^{k,r), (C.5) 
with / = 1,2,3 and a = 1,2,3. Now we can define a further matrix 

n{k,r) = H{k,r)H{k,rbdy)-^ . (C.6) 

This matrix is basically the fields divided by the field's value at the boundary, i.e. if we had 
decoupled differential equations we would get H = diag (A^(r)/A^(rbdy), A^(r)/A^(rbdy), •••) 
Next we take the derivative with respect to r of this matrix, leading us to T-L'[k,r) = 
H'{k,r) ■ H{k,r]^Aj)~^ ■ Again looking at a possible decoupled case we would obtain 
'H'(r) = diag ^A^'(r)/A^ (rbdy), A^'(r)/A^ (rbdy), •• which took the ratio between the 
derivative of the field and the field into account. This is the know result [36,48]. What is 
missing for the Green's function is the r dependent prefactor, which is exactly the content 
of A{k^ r). Therefore we finally get 

G^iu) = 2 lim ^(w, rbdv)^'(w, "Tbdv) + contact- and counter terms, (C.7) 

which allows us to calculate the retarded Green's function in systems where the operators 
may mix. 



^(A;,rbdy) 











bdy 
4 



,B{k,r^Ay) 
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D. General Remctrks on Viscosity in Anisotropic Fluids 



In general, viscosity refers to the dissipation of energy due to any internal motion [49]. 

For an internal motion which describes a general translation or a general rotation, the 
dissipation is zero. Thus the dissipation depends on the gradient of the velocities 
only in the combination = ^ (V^Uiy + S/yU^). and we may define a general dissipation 
function H = \r]^^^^^u^j,yU\p, where r]^^^^P defines the viscosity tensor [13]. Its symmetries 
are given by 



The part of the stress tensor which is dissipative due to viscosity is defined by 

du 



'p.v 



We consider a fluid in the rest frame of the normal fluid = 1. To satisfy the condition of 
the Landau frame U/jJl^" = 0, the stress energy tensor and thus the viscosity has non-zero 
components only in the spatial directions i,j = {x,y,z}. In general only 21 independent 
components of rjijki appear in the expressions above. 

For an isotropic fluid, there are only two independent components which are usually 
parametrized by the shear viscosity t] and the bulk viscosity C- The dissipative part of the 
stress tensor becomes II*-' = —2'q{u^^ — ^6^^u\) — C,u\5^^ which is the well-known result. 

In a transversely isotropic fluid, there are five independent components of the tensor 
jyijfci For concreteness we choose the symmetry axis to be along the x-axis. The non-zero 
components are given by 

^x... = _ 2A , ^yyyy = r?-- = - ^ + , 

^xxyy ^ ^xxzz ^ ^ ^yyzz ^ _ ^ _ ^ (D.3) 

yzyz _ „xyxy _ xzxz _ 

'I — 'lyz 1 'I — 'I — 'Ixy ■ 

The non-zero off-diagonal components of the stress tensor are given by 

n ^ ~ "^VxyUxy ) II = '^fjxy'^xz i ^ 

So far, we only considered the contribution to the stress tensor due to the dissipation 
via viscosity and found the terms in the constitutive equation which contain the velocity of 
the normal fluid u^. In general, also terms depending on the derivative of Nambu-Goldstone 
boson fields = d^ip, on the superfiuid velocity and on the velocity of the director may 
contribute to the dissipative part of the stress tensor. Here the director is given by the 
vector pointing in the preferred direction. However these terms do not contribute to the 
off-diagonal components of the energy-momentum tensor for the following reasons: (1) 
a shear viscosity due to the superfiuid velocity leads to a non-positive divergence of the 
entropy current [49,50] and (2) no rank two tensor can be formed out of degrees of freedom 
of the director if the gradients of the director vanish [43] . In our case the second argument 
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is fulfilled since the condensate is homogeneous and the fluctuations depend only on time. 
These degrees of freedom will generate additional transport coefficients, but they do not 
change the shear viscosities. Thus we can write Kubo formulae which determine the shear 
viscosities in terms of the stress energy correlation functions. 
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